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ABSTRACT 


This thesis investigates the computer aided measurement of fiber 
diameters by laser diffraction. The proposed system consists of a light 
sensitive Random Access Memory (RAM) chip which collects light intensity 
data from the laser diffraction pattern. Measurements of the spatial 
location of the nodes of the diffraction pattern enables the calculation of 
the fiber diameter. These measurements may be performed manually which 
is tedious and requires subjective judgement of the nodes. The alternative 
method of direct processing of the intensity pattern was investigated. 
simulation is conducted to examine the feasibility of this method. Results 
show such a system to be capable of providing one order of magnitude 
greater accuracy than optical microscopy measurements (with a shearing 
eyepiece) and double the accuracy of manual laser diffraction methods 
with the added advantage of permitting the option of total computer 


automation in data interpretation. 
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|. INTRODUCTION 


Fiber reinforced composites are replacing structural metal 
components in today’s aircraft. The high strength and reduced weight of 
composites results in less drag, increased payload, and longer fatigue lite. 
Additionally, the directional properties of composite materials provide 
unique design advantages over conventional materials. 

As with many developing technologies, the reasons for the success of 
fiber composites was not initially appreciated. Tsai [Ref. 1:p. 2] states 
that fortunately the modern composite was so strong it was reliable and 
competitive in spite of less than optimum design practice. Qver the last 
twenty years, much progress has been made in understanding the micro and 
macro mechanics of composite materials. As this Knowledge matures and 
is incorporated into the design process, the full potential of these 
materials can be realized. 

An important contribution to the design process is modeling 
Structural reliability. Development of probabilistic models for a4 
composite must take into account the complex relationships that exist 


between fiber and matrix. 


One important parameter in the probabilistic model is fiber diameter, 
Since variations in fiber diameter affect fiber failure density. Therefore, 
more accurate fiber diameter measurement results in an enhanced 
reliability model and @ more quantitative prediction of structural 
reliability. 

The purpose of this research is to investigate a computer aided 
method of fiber diameter measurement. The existing ‘procedure of 
diameter measurement by laser diffraction is accurate to within 0.5% 
(Ref. 2:p. 210]. It is desired to improve this accuracy by interpreting the 
diffraction pattern with a light sensitive RAM chip. 

The presentation begins with a discussion of diffraction pattern 
analysis. This is followed by an introduction to the Micron Eye's theory 
and operation. Next, the simulation is discussed in depth followed by the 


results and conclusions. 


ll. BACKGROUND 


Fiber diameter measurement by laser diffraction 1s conducted by 
obstructing a collimated laser beam with a fiber sample. A diffraction 
pattern results and is characterized by alternating maxima and minima 


symmetric about a central maximum. (See Figure 1) 
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Figure 1. Fiber Diameter Measurement by Laser Diffraction 


In previous work, Perry, Ineichen, and Eliasson, [Ref. 3] and Bennett, 
[Ref. 4], interpretation of the diffraction pattern consisted of finding the 


distance between interference nodes (minima). This distance is related to 


the fiber diameter. 


Bennett used the slit approximation (Appendix A) to relate the 
distance between nodes to the fiber diameter. Perry, et al., introduced a 
more exact solution by Kerker, [Ref S:p. 260] and compared it to the slit 
approximation. Perry, et al., concluded that the slit approximation should 
be treated with caution. Therefore, Kerker’s solution has been adopted for 
this study. 

In his paper, Bennett successfully demonstrated the feasibility of 
using the light sensitive RAM chip for diameter measurements [Ref 4]. An 
S32 OPTICRAM MICRON EYE, manufactured by Micron Technology, Inc., 
was connected as a peripheral device to an Apple Il+ computer. By 
positioning the Micron Eye so that two interference nodes fit on its 
surface (see Figure 2), an “exposure” of the diffraction pattern could be 
printed by the computer (see Figure 3). Analysis of the printed diffraction 
pattern gave the diameter of the fiber using the slit approximation of 


Fraunhofer diffraction theory (see Appendix A). 





Figure 2. Micron Eye Array with Two Interference Nodes 
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Ht. DIFFRACTION PATTERN ANALYSIS 


Measuring fiber diameters by laser diffraction, requires an 


understanding of the diffraction pattern. This chapter introduces some 


features of diffraction patterns which will be useful in later analysis. 


FUNCTIONAL DESCRIPTION OF THE DIFFRACTION PATTERN 


A. 
The diffraction pattern of a fiber can be described as follows: 


I= (2/KLm) | b,+23b cos(ne) P (1) 


where OQ= the scattering angle 


K, = 20 /A (A = laser wavelength) 


b= Jilad / Ha) 


and a= mdf A (d, = fiber diameter) 


J (a} are Bessel functions of the first kind, 


H (a) are Hankel functions of the second kind. 
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A formal introduction of equation (1) and the related Fraunhofer 
diffraction theory is presented in Appendix A. 

Figure 4 depicts the three dimensional nature of the diffraction 
pattern. Equation (1) describes the intensities along one linear position, 


or slice, of the three dimension pattern. 
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Figure 4. Three Dimensional Diffraction Pattern 


A plot of equation (1) is the Intensity Profile which shows the 
Intensity Ratio (1/1,) versus the angle theta, for a given fiber diameter 


and a given fiber to screen distance L. Figure 5 shows a typical Intensity 


Profile for a fiber diameter of 6.5 microns. 


INTENSITY PROFILE 
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Figure 5. Typical Intensity Profile 
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B. FEATURES OF THE INTENSITY PROFILE CURVES 

The intensity profile has several features worth noting. The first is 
the central (or zeroth order) maximum. The central maximum is the 
largest peak in Figure 35 and it is many times more intense than the 
subsequent maxima. It should also be obvious that the profile is 
symmetric about the central maximum. 

Other features of the curves are shown in Figure 6. Figure 6 is an 
enlarged view of one of the higher order maxima and associated minima. 
As discussed in Appendix A, the higher order maxima are not centrally 
located between the minima, rather they are displaced slightly towards 
the central maximum. Therefore, the higher order peaks are asymmetric 
about their maxima. This asymmetry means the maximum derivative oan 
the upslope side of the curve will be greater than the maximum derivative 
on the downslope side of the curve. This difference is dealt with later 
when tuning for the optimum exposure is discussed. 

The minima in Figure 6 are also distinctive features of the curves. In 
Appendix A, it is shown that the location of the minima are explicitly 


related to the diameter of the fiber using the slit approximation: 


m = * of interference node 
Qi, = sin[ mA/d ] A = laser wavelength (2) 


d = diameter of fiber 
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Figure 6. Features of the Intensity Profile Curve 
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Equation (2) is only an approximation for a fiber, but it is useful for 


preliminary analysis. Such analysis includes: 


(1) determining where to place the Micron Eye for data collection and, 


(2) providing and initial quess of fiber diameter from diffraction 
pattern data. 


C. EFFECT OF DIAMETER ON THE INTENSITY PROFILE CURVES 

Figure 7 shows the effect of diameter on the Intensity Profiles. Note 
that successive maxima and minima are further from the central maximum 
for fibers of smaller diameter. This behavior will be important to Ister 


analysis. 


D. THE MAXIMUM DERIVATIVE AND THE MINIMUM NUMBER OF POINTS 

The goal of automated fiber diameter measurement is to find the 
diameter rapidly and accurately. Speed and accuracy can be exclusive. One 
can imagine that collecting all the points on the Intensity Profile curve 
will result in an almost exact determination of fiber diameter, at the 
expense of time. Two questions must be answered. First, what is the 
minimum number of points to uniquely describe an intensity profile? 


second, where on the curve should these points be obtained? 
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Figure 7. Effect of Diameter on Intensity Profiles (3-D) 


1. Minimum Number of Points 

Equation (1) consists of an infinite series. Determining the 
minimum number of points to uniquely define an infinite series is not 4 
Straightforward. procedure. Cursory inspection of Figure & leads to the 
conclusion that one point will not uniquely define the curve, since 
intensity curves for many diameters pass through the same point (e.g., at @ 
= O4radians). If the absolute intensity is not known, 2 points will not 
provide 4 unique solution either. For the case of non-absolute intensity 
measurements (as is the case in this experimental set-up) 3 points 
appears to uniquely define a curve. 

The curve in Figure 9 illustrates this observation. An imaginary 
line is fixed through three points (i.e., the experimental measurements). 
The absolute intensities of these three points are not known, and neither 
- are the absolute spatial locations with respect to the central maximum 
(ie, ©). Only the relative spatial locations among these points are 


known. Iteration of the diameter is equivalent to moving this line (with 
the x marks fixed relative to each other) along different locations on the 
curve, trying to match all three points. This iteration can be repeated on 
curves of different diameters (Figures 10 and 11) and no other match can 
be found. Thus, this argument defines three points as the minimum 


required to uniquely determine the intensity profile curve. 
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Figures 9, 10, 11. Three Points to Define a Curve 
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2. Maximum derivative 

Having heuristically established that three points are required, the 
next step is to determine the optimum location from which to select the 
points. For example, if points are selected at the mimima for maxima) a 
wide variation in theta results from a small change in Intensity Ratio. 
This variation is defined by the derivative of the Intensity Ratio with 
respect to Theta, and at the extrema this derivative is very small. It can 
be shown that the points on the curve where the derivative is a maximum 
will have the least variation error. Therefore, it is desirable to use the 
Intensity Ratios corresponding to the maximum derivatives as the 
optimum intensity ratios for the three points. These Intensity Ratios will 
be referred to as the Threshold Intensity Ratios. 

Figure 12 shows a portion of an Intensity Profile with a derivative 
curve (absolute value) superimposed. ‘Yertical lines drawn through the 
derivative curve maxima intersect the Intensity Profile showing the 


location of the optimum (threshold) intensity points. 
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Figure 12. Intensity Profile with Superimposed Derivative Curve 
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IV. MICRON EVE 


The MICRON EYE image sensor is an optically sensitive Random Access 
Memory (RAM) chip capable of sensing an image and translating it to 
digital computer compatible signals. The Micron Eye was selected 
introduce automation to the existing techniques of fiber diameter 
measurement by laser diffraction. Automation is expected to increase the 
speed and accuracy of diameter measurements. This chapter introduces 


the theory and operation of the Micron Eye. 


A. PHYSICAL LAYOUT AND DIMENSIONS 
The Micron Eye (1S32 OpticRAM ™) has two arrays each containing 126 
rows x 256 columns of sensors. This application will use only one of the 


arrays. The size of an array is 4420 ym by 877m. (See Figure 13} 
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Figure 13. Micron Eye Array 
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Each sensor is a light sensitive element is called @ pixel. The 
physical organization of the pixels is shown in Figure 14. The 126 x 256 
elements actually map into a 129 x 514 “cell placement grid”. This 
arrangement leaves “space pixels” in between each pixel in the row 
direction. The space pixels can be set high, low, or to the level which 


agrees with the majority of its nearest neighbors. 


B. THEORY OF OPERATION 

The pixels are capacitors which discharge 4 preapplied voltage at 4 
rate proportional to both the intensity and duration of the impinging light. 
The voltage in an exposed capacitor is read and digitally compared to the 
fixed threshold voltage. If the voltage is below threshold the pixel is read 
as WHITE. If the voltage is above threshold the pixel is read as BLACK. 
The digital comparison concept will be an important part of the simulation 
in the next chapter. 

After a pixel is read, the row containing that pixel is refreshed. 
Refreshing sets all pixels which are below threshold to 0 volts, and all 


pixels which are above threshold to +5 volts. ([Ref. 6] 
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Figure 14. Micron Eye Physical Organization 
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1. Operation of the Micron Eve 
Operation of the Micron Eye is simple. The current research 


configuration uses an Apple Macintosh (512K) computer. Control software 
is provided by Micron Technology, Inc. The Micron Eye is connected 
through either the modem or printer port. (See Figure 15) The computer 


also acts as the power source for the Micron Eye. 





Figure 15. Macintosh and Micron Eye 


2. Exposure 
An exposure of a portion of the diffraction pattern can be made by 
varying exposure times. A sample exposure is shown in Figure 16. It is 
Important to visualize that this exposure represents across section of the 


intensity profile curve. 


intensity Profile 


Micron Ege 


Intensity Ratio 





Figure 16. Exposure and Its Relation to Intensity Profile 


Longer and shorter exposure times will vary the size of the cross 
section. Yarying exposure time is equivalent to moving up or down the 


intensity ratio axis of the intensity profile curve. 


C. INTENSITY MEASUREMENTS WITH THE MICRON EYE 

Bennett [Ref. 4] used the Micron Eye to measure the distance between 
interference nodes. This approach was compatible with the slit 
approximation requiring only the theta locations of the interference nodes 
to give the diameter. In Kerker’s equation [Ref. S:p. 260], one cannot solve 


explicitly ford. The solution requires Knowledge of the Intensity Ratios 


2 


and their respective theta locations. For the Micron Eye, Intensity Ratios 
correspond to the user controlled exposure time. 


1. Two Approaches to Intensity Calibration 


The Micron Eye must be calibrated to a reference intensity level. 
One method to accomplish this would require two exposures. One at I, + 
A and another at I, + A. Delta represents an unknown level above the 


zero (or absolute) intensity. (See Figure 17) The difference between 
these two values eliminates delta. 
A second method requires a fiber of known diameter. Qne exposure of 


this “calibration” fiber will allow determination of the absolute intensity. 


INTENSITY PROFILE 


Figure 17. Calibrating the Micron Eye 
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¥. SIMULATION 


This chapter introduces the simulation of the Micron Eye / Laser 
Diffraction diameter measurement system. The simulation answers many 
questions relating to the real system: 

(1) where should the Micron Eye be placed? 
(2) what exposure is best? 


(3) what accuracy can be expected? 
(4) is the computer code valid? 


The simulation combines diffraction pattern analysis with Micron Eye 
operation. The worth of the results will depend on the accuracy of the 


simulation. 


A. COMPUTER PROGRAMS 
several computer programs have been written to conduct this 


Simulation. These are: 


(1) DATAMAKR 
(2) EXPOSURE 
(3) DIAFIND 


Appendix B discusses these programs in detail. 
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Briefly: 
(1) DATAMAKR is used to produce Intensity Ratio versus Theta data 
for a given fiber diameter d, and screen to fiber distance L. The 
user controls the number of points produce and their spacing. 


(2) EXPOSURE reads the Intensity Profile data generated by 
DATAMAKR and recommends the optimum exposure based on the 
average of the intensity ratios corresponding to the maximum 
derivatives. Exposure will also introduce random error into the 
data, as desired. 


(3) DIAFIND uses the data generated by EXPOSURE to find the diameter 
of the fiber. The programs begins with 4 user input guess of the 
diameter and conducts an iterative search / comparison routine to 
find the actual diameter of the data. 


B. OVERVIEW 

The general approach in examining the proposed system will be to consider 
a typical carbon fiber. This fiber has been assigned the arbitrary diameter 
of 7.254u1m. Two additional fibers +20% of the 7.254um fiber are also 
considered. These fibers define the range of diameters from 3.603um to 
G./O3ym. See Figures 16, 19, and 20 for Intensity Profiles curves for 


these diameters. 
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Figure 18. Simulation Intensity Profile (5.603 1m) 
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Figure 19. Simulation Intensity Profile (7.2541m) 
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SIMULATION INTENSITY PROFILE 
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Figure 20. Simulation Intensity Profile (8.705u1m) 


C. POSITIONING OF THE MICRON EYE 

In order to measure the diameter of the typical fiber, a method of 
postitioning the Micron Eye must be defined. Further, from this same 
position, it is desired to also measure the +20% fibers. This requires that 
one position of the Eye permit fiber diameter measurements over the 


specified range of fiber diameters. 


oy 


Such @ position is defined such that three data points will be obtained 
for any fiber in the diameter range, as in Figure 21. First, a simpler case 


will be discussed to introduce positioning for a given diameter fiber. 


intensity Profile 


Intensity Ratio 





1 2 3 Theta (radians) 





Micron Eye 


Figure 21. Exposure for Three Data Points 
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1. Pasitioning for a Single Fiber 


suppose an exposure of the first through the second interference 


nodes (2°! maximum) was desired for the 7.254ym fiber. Figure 22 shows 


the relationship between the fiber and the Micron Eye: 


FIBER 





Figure 22. Micron Eye Placement in Diffraction Pattern 
@. = tam’ [x /L] (3) 
6, =tan'[x,,/L] (4) 
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or tan Giti-tan Gi = U/L (xitt - xi) = Ax/L (5) 


take Ax = 4.4 mm (width of Micron Eye array) 


With Ax fixed, one can adjust L so that any desired O.., - 0. will 


Recall that @.. = sin [ mA/ d ] 


fit on the Micron Eye's array. — 


equation (2), which will give the theta locations of the first and second 


interference nodes. 


Onin = Sin (mA/d) | where m=node * 


A = laser wavelength (632.8nm) 


d = fiber diameter 


for 7.254y1m: 0. = 086712 radians 


O., .17459 radians (fori=- 1} 


with ©. and 0.,, fixed, L is defined from equation (5): 


L =- ax # [ tan 6.,, - tan 0. | = 49mm 
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L is the distance of the Micron Eye from the fiber in the 
longitudinal direction. Referring back to equation (3), the lateral 
placement of the Eye is given (see Figure 23): 


x - L tan{e@) = 427 mm 


inner 


FIBER 





Figure 23. Micron Eye Position 


Thus, equations (3) through (5) give the Micron Eye position for any 
desired d, 0. 


1#+1? 


0., L and m. 
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This position is not yet optimum for single diameter case. This 
analysis has captured only two interference nodes and the intervening 
maximum. Any exposure will yield information somewhere between the 
two extremes shown in Figure 24. The Micron Eye images correspond to 
the shaded areas of the intensity profile curve on the right side of Figure 
24. 
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Figure 24. Short and Long Exposures 


One method of finding the points would be to search the Micron Eye 
array data for a major axis. The endpoints of this major axis are the 


points required for the analysis. (See Figure 25) 
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Major Axis 





endpoints 


Figure 25. Endpoints of Major Axis 


Recall from diffraction analysis (Chapter 3) that three points are 
required to define the intensity profile curve. Therefore, more of the 
diffraction pattern must De seen by the Micron Eye. A reductionin L will 
yield the third point. (See Figure 26) L should not be reduced any more 
than is necessary, since fewer pixels are being used to describe the data, 


degrading the resolution of the Micron Eye. 


Major Axis 





endpoints 


Figure 26. Three Points on Major Axis 
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2. Positioning the Micron Eye for _a Range of Diameters 
Consider the positioning of the Micron Eye so that three data paints 


can be collected for a range of diameters. A “window’ for the array must 
be defined which will ensure three points of data for any diameter in the 
specified range. 

Examine the Relative Intensity Profile plot in Figure 27. © begins 
at .0735 radians which excludes the the central maximum. The central 
maximum is so intense, the Micron Eye’s fastest exposure cannnot prevent 
overexposure. Therefore, the central maximum is not considered as 4 
location for the window. 

The window should also be located in 4 region where the maximum 
intensity corresponding the the largest diameter fiber is close to the 
maximum intensity of the smallest diameter fiber. This ensures the 
Threshold Intensity Ratio will be more representative for all diameter in 
the range. This condition occurs at the higher order nodes (2 or greater), 
as illustrated in Figure 26. The plot of intensity ratio derivatives also 


shows this condition in Figure 29. 
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Figure 27. Relative Intensity Profiles 
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Figure 28. Region of Similar Maximum Intensities 
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Figure 29. Region of Similar Intensity Ratio Derivatives 
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To obtain three points, at about 1.5 maxima are required. Since 
the distance between nodes is the greatest for the smallest diameter, find 


the @ locations for the smallest diameter which gives 1.5 maxima. 


for 5.803 um: 
Owing = SiN (21/5.803pm) = .216 rads 
0 ina = Sin (31/5.803pm) = .321 rads 


Half of this range is about .050 radians and 1.5 maxima can be 
approximated by .165 ----- >» .321 radians. Admittedly, this is not 4 
sophisticated method. The method is somewhat liberal and it was later 
discovered that a range of .165 to .300 radians was better than the larger 
range. This is because the minima are not of interest so the right hand 
minimum at .321 radians was discarded, and the range reduced to 300 
radians. Only the part of the curve where the slope is a maximum needs to 
be “seen” by the array. 

A quick check can be made to insure that the curve representing 
the maximum diameter in the range fits inside the window. 


for 8.705pum: 
Ono = Sin (21 / 8.705) = .150 rads 


0 nin’ = sin (3 1/ 8.705} = 216 rads 


@_ 4 = Sin (41 / 8.705) = 286 rads 


4/7 


The next step is to determine the distance L at which the Micron 
Eye will cover the @ range of .165 


---> 300 radians. 
defines the 


longitudinal position and equation 
corresponding lateral position: 


Equation (5) 


(3) defines the 


geass 
I 


Ax? [tan(6., .) - tan(9.)] = 30 mm 


on 
i 


, =L tan(o.) « Smm (inboard position of array) 
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Vi. SIMULATION: THE PERFECT AND IMPERFECT DATA SETS 


Now that a window has been defined (0 range) for the diameter range, 
it is necessary to construct perfect data tor the simulation diameters: d,, 
d,, andd,. DATAMAKR generates this data for 4 given d and L, over the 
specified range of 9. Additionally, the number of points over the specified 


range must be chosen. The number of points defines the © interval and 


corresponds to the physical spacing of the pixels in the Micron Eye array 
which 1s on the order of 104m. To match this spacing, 400 points were 


generated over the interval .165 to .300 radians. 


A. DIGITIZATION OF DATA 

The EXPOSURE program is used to digitize the perfect data. Digitizing 
the data is the simulative analog of the Micron Eye threshold voltage 
comparison. Inthe simulation, any intensity ratios above a certain ratio 
(call it the 7Aresha/d retja) are assigned a value of 1 and any intensity 


ratio below the threshold ratio are assigned a value of 0. 


49 


1. Determining the Threshald Intensity Ratio 
Before EXPOSURE can digitize the perfect data, the threshold ratio 


must be determined. EXPOSURE calls the subroutine DERIV which 
calculates the derivatives of the intensity profile data at each theta 
location. 

DERIVY proceeds to search for all the local maximum derivatives. 
The average of the maximum derivatives is returned to the EXPOSURE 
program and is used as the threshold intensity. 

After the data is digitized with respect to the threshold intensity, 
the digital data is searched to find the theta locations where the intensity 
changes from 0 --> 1 or 1 --> 0. The values for theta at these locations 
are averaged and this average is taken to be the theta location where the 
threshold intensity is located. The averaging of theta is based on the fact 
the Intensity Profile curve is approximately linear in the region of the 
Threshold Intensity Ratio. Figure 30 shows the general location of the 
Threshold Intensity Ratio. 

In the simulation, EXPOSURE finds three average theta locations 
corresponding to the threshold intensity ratio. These data are icred in 4 
data file which is subsequently input to the program DIAFIND. DIAFIND 
prompts the user to input an initial guess of the diameter and proceeds ta 
find the diameter associated with the EXPOSURE data. 
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Figure 30. Averaging of Theta 
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B. INTRODUCTION OF ERROR INTO THE PERFECT DATA 
EXPOSURE prompts the user for an error input. If error is desired, 
EXPOSURE calls the subroutine RANDOM. RANDOM introduces error based on 


the maximum intensity ratio in the perfect data: 


l-=lp + I I, * ERROR * RND | 
where ERROR is the error to be introduced and RND is a random number 
between +1 generated by the NONIMSL subroutine RANDU. RANDOM returns 


the now “imperfect” intensity ratios to the EXPOSURE program. 


C. DIGITIZATION PROBLEMS WITH THE IMPERFECT DATA 

Digitizing the imperfect data results in local regions where the 
intensity oscillates between O and 1, as shown in Table 1. The 
introduction of random error changes the smooth curve to erratic points. 
Locally, variations above and below the threshold intensity ratio cause a 
series of O-->1 and 1-->0 oscillations. (See Figure 31) Ideally, only one 
local theta location is to be associated with the threshold intensity ratio. 

There are two approaches to this problem. First, the average of the 
local theta locations can be calculated. This method results in three 
theta locations for input into the DIAFIND program (the same as perfect 
data). 


az 
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TABLE I, 


. 2321624000E 
. 2324997000E 
. 2328373000E 
. 2331748000E 
. 2335124000E 
. 2338498000E 
. 2341873000E 
. 2345249000E 
. 2348623000E 
. 2351998000E 
. 2355374000E 
. 2358747000E 
. 2362123000E 
. 2365499000E 
. 2368872000E 
. 2372248000E 
. 2375622000E 
. 2378997000E 
. 2382373000E 
. 2385748000E 
. 2389124000E 
. 2392498000E 
. 2395873000E 
. 2399249000E 
. 2402623000E 
. 2405998000E 
. 2409374000E 
. 2412747000E 
. 24161230005 
. 2419497000E 


DIGITIZATION ERROR 


THETA LOCATION 


INTENSITY RATIO 


=> Gao CO CO @,e2eo2> Cie @ © GeO © OO @ Oe © 2 OO 0 @ 0 CO Ge Oo 


- 3544308000E-04 
- 3196243000E-04 
. 3160309000E-04 
- 3246925000E-04 
. 3372986000E-04 
. 3785736000E-04 
. 3090432000E-04 
. 3342219000E-04 
-3387519000E-04 
- 92053730005 -04 
- 35 162220005-04 
. 3860490000E-04 
- 3421140000E-04 
. 3468163000E-04 
- 3580747000E-04 
- 3904779000E-04 
9) 7 G2 7OUUUE=04 
» 3912353000E-04 
. 3599435000E-04 
- 3474336000E-04 
- 3538676000E-04 
- 3470230000E-04 
- 3952176000E-04 
po J2o007 VOOR 04 
-3944175000E-04 
. 3966210000E-04 
-4409726000E-04 
-4269459000E-04 
-4371968000E-04 
. 3904802000E-04 
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Figure 31. Intensity Variation at 5% Error 
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The second method is to accept the theta locations as they are and 
enter them into the DIAFIND program. The latter method will pass a 
variable number of data points to DIAFIND. 

Both approaches were tested (at 5% error for 7.254 jum) and the same 
diameter was recovered by DIAFIND in each case. The only difference 
between the two methods is that the averaging method runs one second 
faster (out of 22 seconds on the IBM 360) than the other method. 

The averaging method was adopted for this simulation for two 


reasons: 
(1) The data input to DIAFIND will always be a constant number of 
points so that DIAFIND can be easily adapted to run actual Micron 
Eye data without simulation related logic buried in the code. 


(2) The averaging method returns the same diameter as the other 
method. 


D. CHOOSING THE CORRECT EXPOSURE (TUNING) 

The optimum exposure for this simulation was calculated by the 
program EXPOSURE. It was an average value of the intensity ratios 
associated with the maximum derivatives locations for the ?.254um fiber 
data. This same exposure was used throughout the simulation for all the 
fibers. 

Using the same exposure introduces realism into the simulation since 
one cannot tune the exposure for every fiber that is to be measured. 


Because the Micron Eye window is located at least two nodes out from the 
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central maximum the exposure is closer to optimum for all diameters over 
the range. {In contrast to a window location closer to the central 
maximum.) 

The process EXPOSURE uses to determine the threshold intensity ratia 
can be called “tuning”. The Micron Eye analog of tuning would consist of 


three steps: 


(1) Determine the approximate diameter for the fiber to be measured 
using an optical shearing eyepiece. 


(2) Determine the Micron Eye location parameters L and x by the 
methods in Chapter 5. 


(3) Vary exposure to obtain a defined relationship between the three 
points on the Micron Eye photograph. 


AS an example, return to the case where L = 30 mm and Kener = smm. 


(x. _., is the inboard x location of the Micron Eye array.) Figure 32 shows 
the “tuned” relationship between three points on the Micron Eye array. 
This relationship can be expressed as a ratio of a:b:c. The theta 
locations at the inboard and outboard edges of the Micron Eye are known. 


e.g., One = Oo radians and 6 = .300 radians. 


outer 
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0, @ Oo, 0 outer 


Figure 32. Three Point Spacing Ratio 


This means that 2 Pea 0, < 0, « Oz < oe 


To find 6, , ©, , and Oz. run EXPOSURE for the diameter being tuned. 


EXPOSURE outputs the theta values for the optimum exposure. (See Figure 
33) 
The theta values are related to x locations by: 


x, =L tan(9, ) 
Xo =L tan(o,,) 


Re =L tan(@,) 
where a=, -k 

= Ro ~Xy 

E> tails 


of 


intensity Profile 


ir 





8 (radians) 


Figure 33. Finding the Theta Locations 
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INTENSITY YR AITO = 10 


0.03.0 6.0 9.0 12.0 15.0 18.0 


YIl. DISCUSSION OF RESULTS 


The simulation was conducted for three diameters: 


d, = 9.603 yim 
d,= 7.254 jim 
d, = 8.705 um (where d, and d, are +202 d.) 


For each d,, three levels of error were introduced: 1%, 2% and 52%. 


Figure 34 shows the three intensity profile curves for the d.. 


SIMULATION INTENSITY PROFILES 


LEGEND 
7.254 MICRONS 
9.803 MICRONS 
8.705 MICRONS 





OL075 O20 0.175 O2225 Ocevo Qeaco Clyeslve 
THETA (RADIANS) 


Figure 34. Simulation Intensity Profiles 
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The simulated Micron Eye window was determined as previously 
described at L = 30mm and ner = anim. The exposure was tuned with 


respect to the ?7.254um perfect data. EXPOSURE recommended the 
Threshold Intensity Ratio of 0000370. This Threshold Intensity Ratio 
remained constant throughout the simulation, for all diameters. 


Before the simulation, all diameters were tested with no error using 
the 7.254ym Threshold Intensity Ratio. DIAFIND recovered d,, d, and d, 


exactly (i.e., to the three digits accuracy of the original diameters). 

The simulation originally began by collecting thirty data points for 
each diameter/error combination. Twenty additional points were 
collected (50 total) to produce meaningful histograms. 

Figures 35 through 423 are histograms depicting the results of the 
simulation. In general, the expectation was to see a decrease in 
resolution as more error was introduced. It was also anticipated the 
resolution would be best for the 7.254um case (for all values of error) 
since the exposure was “tuned” for this diameter. Finally, it was hoped the 
method would be more accurate than existing laser diffraction 


measurement methods (= 0.53%). 
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Figure 35. 5.803pm Fiber (12% error) 
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7.254 fim fiber (1% error) 
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Figure 36. 7.254ym Fiber (1% error) 
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8.705 ym fiber (1% error) 
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Figure 37. 8.705y1m Fiber (1% error} 
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Figure 38. 5.803y1m Fiber (2% error) 
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Frequency 


7.254 pm fiber (2% error) 
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Figure 39. 7.254ym Fiber (2% error) 
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6.705 pm fiber (2% error) 
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Figure 40. 8.705ym Fiber (2% error) 
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Frequency 


5.803 jm fiber (5% error) 
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Figure 41. 5.603y1m Fiber (5% error) 
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Figure 42. 7.254j1m Fiber (5% error) 
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Figure 43. 8.705jm Fiber (5% error) 
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A. ACCURACY VERSUS RESOLUTION 

The histrograms show the accuracy and the resolution of the method. 
The accuracy is associated with the largest spike, and related to the size 
of the interval called the resolution. The resolution is the ability of the 
routine to “see” a difference between two different fibers. For example, if 
the resolution is .0O5um the method does not discriminate between 


6.705y1m and 6.700um. Table 2 shows accuracy versus resolution for the 


results. 
Table 2. ACCURACY VERSUS RESOLUTION 
SERROR DIAMETEREM RESOLUTION ACCURACY (Res/dact) 2 
1 5.603 004 069 
1 2a 0042 056 
6.705 9060 069 
Z 5.603 0046 083 
Pp fay. | 0054 O74 
6.705 0049 056 
5 5.603 0105 161 
Eea4 0070 096 
6.705 0136 .156 


70 


The data show a decrease in accuracy and resolution as error is 
increased. It is evident in most cases that the Tesco results are 
better because it was the “tuned” diameter. The only exception is the 
accuracy for the 8.705um fiber (2% error) is better than the 7.254ym 
fiber. Overall, the largest error is .18 percent which is less than one half 


of the error associated with the manual laser diffraction method. 


B. ERROR 

There are two contributions to error in this simulation. 

The first can result from programming/calculation errors. Many trial 
runs of the software were taken to minimize the likelihood of this kind of 
error. 

The next type of error is the digitizing error encountered in an actual 
experiment. This error is simulated based on the maximum intensity of 
the profile curve, in the region of interest, to apply random error equally 
to all points. Had the random error been based on each point, points with 
less intensity would have less error, and points with higher intensity more 
error, which is inconsistent with the physical digitizing process. Using 
the maximum intensity avoided such a condition but it cannot be 


ascertained that this is the optimal representation of the physical system. 
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Yill. CONCLUSIONS 


The results of the simulation are encouraging. The largest error is 
18% or less than one half of the manual methods. Because the programs 
were carefully developed and tested it is unlikely they contributed to the 
error. Also, much effort was directed towards accurate simulation of the 
physical system so that the results would reflect what can be expected 
from the actual experimental measurements. 

The simulation demonstrates an increase in accuracy two to ten 
times better than that currently possible by making manual measurements 
with laser diffraction. The method also lends itself to automation which 
makes it attractive for quality control purposes and research for 


materials development. 
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IX. RECOMMENDATIONS 


This study has shown the feasibilty of computer aided diameter 


measurements. There are many directions future work can take. 


4. SOFTWARE. 

More development and testing of the software could result in 
increased accuracy. It would also be valuable to implement the software 
on a small computer (like the MacIntosh) to allow real time processing of 


actual data. 


B. HARDWARE. | 

There remain some hardware considerations which must be resolved. 
The biggest of these, perhaps, is the accurate positioning of the Micron Eye 
array in the direction perpendicular to the laser beam. In an effort to 
increase the resolution and accuracy of the problem, the system may 


benefit from two Micron Eyes. 
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APPENDIA A. FRAUNHOFER DIFFRACTION THEORY 


The following is 4 brief discussion of Fraunhofer diffraction theory 


with emphasis on aspects of the theory which relate to this research. 


4. FRAUNHOFER DIFFRACTION | 

Fraunhofer diffraction (Figure 44) results when light approaches and 
leaves an obstacle or aperture in the form of plane wavefronts. [Ref. 4:p. 
176] The light source and the plane of observation in effect are at 
infinity. A collimated laser beam is ideally capable of Fraunhofer 


diffraction because its beam consists of parallel rays advancing in phase. 


Central Maximum Interference Nodes 


wi 





Figure 44. Typical Single Slit Diffraction Pattern 
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B. THE CLASSICAL SINGLE SLIT EXPERIMENT 

The simplest demonstration of Fraunhofer diffraction is the single 
slit experiment. (See Figure 45) Parallel, collimated light passes through 
aslit of width a. The diffraction pattern is visible on a screen located 6 
distance L from the slit. An observer at point P, moving across the 
Seneen, sees & succession of maximum and minimum intensity points. 
These extrema are the result of constructive and destructive interference 
of the light. For example, a minimum occurs when the angle © produces 4 
phase difference of one wavelength between the rays at the upper and 


lower edges of the slit. Thus, minima occur whenever 
asin(O@)=mA, wherem= 1, 2,3..... (A.1) 


These minima are referred to as interference nodes. (For further 


discussion of this subject, see Meyer-Arendt, [Ref. 4].) 
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Figure 45. Single Slit Path Difference Relation 
to the Interference Node 


1. Diffraction Minima 
Examination of equation (A.1) shows that as the slit becomes 
narrower the angle © becomes larger. As this theory is extended to 
approximate the diffraction phenomena of an obstruction (a fiber), one can 
expect short distances between interference nodes for larger fibers and 


greater distances between nodes for fibers of smaller diameters. 
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Another important point concerns the distance between minima. It 
appears the interference nodes are equidistant. This is true only within 
the limits of the small angle approximation. The distances between 
Interference nodes actually increases as one moves outward from the 


central maximum by the relation: 


OG uig = Sin [ mA/d | (A.2) 


2. Diffraction Maxima 
The diffraction maxima are not located midway between minima. 
The locations of the maxima can be derived as by Meyer-Arendt [Ref. 7?:p. 
220]. These occur whenever the derivative of the intensity is equal to 


Zero: 


di/dB = di/dB [1, (sin B/B8)7] = 0 


2lo sin8/8 [ -sinB/B* + cosB/B | 


or whenever: B = tanB (A.3) 


t? 


As per Houstoun [Ref. 8:p. 68] the maxima can be seen where 
y=tanB and y=B. (See Figure 46) The maxima are displaced slightly 
from center towards the central maximum. Much further away from the 


central maximum, the maxima ere nearly halfway between minima. 


Figure 46. Locations of the Diffraction Maxima 
C. FRAUNHOFER DIFFRACTION AND THE FOURIER TRANSFORM 


The foregoing discussion centered on the Fraunhofer diffraction due 
to a slit. The mathematics of the slit example are simple and provide 
insight into diffraction physics. A more elegant derivation of Fraunhofer 
diffraction shows that the diffraction pattern of an object is the Fourier 


transform of that object. [Ref. 9:p. 174] 
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for the slit: 


fx) 









| 0, Ixl>b 
f(x) = Be eh ek 


F(O) = 2 sinbo 
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Figure 47. Fourier Transform for the Single Slit 


Thus, F(@) describes the amplitude of the diffraction pattern, and 
[F(@)]* represents the intensity. 


Modeling the transform for the obstruction is more complicated. 


consider the complement of the slit transform: g(x) = 1 - f(x). 


g(x) = 1-f(x) _ Ft, Ixl>b 
g(x) = bs Ix] < b 


SSS ee EL IE) 





-b +b a 


Figure 48. Complement of the Fourier Slit Transform 
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This results in the form [ oe - 2sin (b@}/O ] which is not 
transformable. | 
An alternative is to use a substitute function h(x) = g(x) - f(x). 


Physically, this is approximating the obstruction as two parallel slits: 


h(x) = g(x)-f(x) 





-C -b +h +C 


Figure 49. Fourier Transform for Two Parallel Slits 


which has the solution: 
H(@) = (2sin c90)/O - (2sin bO)/ OC 


where c would be large, but not infinite. 


This is still a crude approximation which shall be improved upon in 


the next section. 
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D. REJECTION OF THE SLIT APPROXIMATION 

In their paper titled “Fiber Diameter Measurement by Laser 
Diffraction” [Ref. 3:p. 1378], Perry, Ineichen, and Eliasson conclude that 
the diffraction pattern of a real fiber is sufficiently different from that 
of a slit to warrant the slit approximation being treated with caution. 

Further, they recommend 4 solution presented by Kerker [Ref.5:p. 260] 
which has been adopted in this study. Kerker's solution assumes the fiber 
is perfectly reflecting (i.e.,the reflective index m=co), and while this is 
not completely true, Perry, et al., [Ref. 3:p. 1378] indicate some degree of 


absorption is not likely to be significant. 


E. INTENSITY EQUATION FOR A REAL FIBER 
Kerker [Ref.5:p. 260] gives the scattered intensity relation for a real 


fiber: 


i= (2/K Lx) | b,+23b cos(ne) F (A.4) 


where O= the scattering angle 


K, = 20 /A (A = laser wavelength) 


b= Ja) / Ha) 
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and a= md,./1 idee fiber diameter) 
J (a) are Bessel functions of the first kind, 


H_ (a) are Hankel functions of the second kind. 


The real fiber equation (A.4) is somewhat obscure in its compact 
form. It can be shown that: 


1/1 =(2/K Lo) [« r.+22 r.cos(nd) )7+(s +25 s_cas(nd) 2] (A.5) 


where r 


Jo) J Ld2a) + ¥ 2a) |} (A.6) 


and s, = Ja) ¥ lo) / [db 2ad+¥ 2a] (A.7) 


Now, the intensity ratios for any @ location can be calculated for any 


diameter fiber. 


l. Sensitivity of the Results to the Number of Bessel Terms 
The calculation of Intensity Ratios requires the computation of n 


r, and s| Bessel terms, equations (4.6) and (A.7). Here the phrase 


“Bessel terms” indicates the algebraic combinations of the J, and ¥, 


| Bessel functions. 
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a. Two competing phenomena 

There exist two competing phenomena which govern the number 
of Bessel terms to be used in the calculations. The first requires a 
minimum number of terms for accuracy. The second, limits the number of 
terms so the Y functions do not cause an underflow error during 
computation. 

(1) Minimum Number of Terms. As in any series, there is 4 
minimum number of terms required for computational accuracy. Figure 30 
shows the effect of the number of Bessel terms computed for seven 
curves, all with a diameter of eight microns. 

Note that the 43, 50, 75, and 86 term curves are nearly identical and 
the curves with fewer than 43 terms are decreasing towards a smooth 
line. This research, indicates that 50 terms returns four digit accurecy 
for diameters ranging from 5 to 10 microns. 

(2) Maximum Number of Terms. The maximum number of terms 


depends onthe value ofa. As the diameter decreases,so dcesa@. This 
in turn produces very large values of ¥ 4{a) (in the fF and s. 


denominator} which causes an underflow error. For the diameter range 
considered in this research, underflow was not a problem. In one 
instance, below 5 microns, it was found the number of Bessel terms had to 
be reduced tc 43 to prevent underflow. The effects of the minimum 


number of terms was not investigated at this smaller diameter range. 
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Figure 50. Effect of Number of Bessel Terms 


2. Effect of Different Diameters 
Recall that the diffraction pattern for a slit showed that as the 
diameter decreased, the interference nodes moved away from the center of 
the pattern. It is interesting to note that the pattern for a real fiber 


behaves the same way. (See Figure 51) 
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Figure 31. Effect of Different Diameters 


APPENDIX B. COMPUTER PROGRAMS 
The following is a list of programs discussed in this Appendix: 


DATAMAKR 
DIAFIND 
EXPOSURE 


The programs are written in Waterloo Fortran I'v (WATFI'Y) and run on 


the Naval Postgraduate School's IBM 360 computer. 
4&4. FORMAT OF DATA 


All data is formatted using exponential notation. The most 
frequently manipulated files are those containing theta locations and 
intensity ratios. The format for these files is: (1%, £17.10, 1X, £17.10). 
In WATFIY all data files must be of filetupe “WATFI¥”. 

To compile and run a program on the IBM 360 type “© WATFI¥ 
PROGNAME DATAFILE *{XTYPE", where PROGNAME is the filename of the 
program and DATAFILE is the filename of the data file. 
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B. COMPUTER PROGRAMS 
This appendix presents the various computer programs and outlines 


their logic. 


1. DATAMAKR 
DATAMAKR generates the data used in the simulation. The equation 


for computing the Intensity Ratios of the diffraction pattern for any value 


of theta is: 
l= (2/K Lm) | b,+22b, cos(né) ‘a (A.4) 
where 6= the scattering angle 
9 = 20 fA (A = laser wavelength) 
b= Ja) / Hla) 
and a= Wd. /A (d, = fiber diameter) 


J (a) are Bessel functions of the first kind, 


H 2a) are Hankel functions of the second kind. 


All computations Begin by calculating the required number of 


Bessel terms. Two HONIMSL subroutines are called: BESJ for the J ba) and 
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BESY for the ¥ fa). Note: should double precision be desired, the NONINSL 


subroutine BINT will return double precision values for the J. and ¥ . 


The Bessel function values are stored in an vector JI) and ¥(1), 
and are used to calculate values of R{1) and S(I). Next, the matrix of 
cosines is contructed. The size of this matrix is M by K, where N = M-1 
and N is the number of Bessel terms. K is the number of theta locations. 
Values of theta at each location have previously been stored in an array 
called T(K). 


To visualize the program's calculations, the matrices symbolic of 


these calculations are sketched. (See Figure 52) Note that r, ands, are 


outside the summation term in equation (4.1) and therefore vill not be 
multiplied with any cosine terms. Since array subscripts must be 


designated with a nonzero value, rr, ands, will be represented by r{1) 


and s¢1). 

After multiplying the two matrices, it is necessary to complete 
the summation by summing the columns in the matrix. This results in an WN 
term value for each value of theta. A straightforward calculation then 
yields the intensity ratios. The intensity ratios are stored in an array 


with their corresponding theta locations. 
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r, coso, cos 6,5 ee 
rx cos20, cos 20,5 
a x | £8836; cos 30. 
re cos(m-1)6, cas(m-1)6, 
La a a 
Fa" 
[R] x {C] _ rc. 
tn om-1 ly Ae si 
im 
* r,cos(n6 ) 
n=1 


= = * 


tn“ m-1 Kk 


Figure 52. DATAMAKR Matrix Algebra 
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The logic of DATAMAKR is: 


£9] 8)0 tetas > * Bessel terms to be computed 
diameter of the fiber 
screen to fiber distance, L 


array of theta location values 


Calculate: R, and S using BESJ and BESY 


Produce the Matrix of Cosines 
Produce the Matrix of Products 
Sum the columns of the Product Matrix 


Compute and output the Intensity Ratios 


2. DIAFIND 





This program finds the diameter of a fiber through an iterative 


process of residual comparison. The program accepts the data which has 


been output by the program EXPOSURE. The user is prompted for K and d.. 


K is the number of theta values the program should expect and d. is the 


initial guess of the fiber diameter Cin microns). 
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The program first calculates intensity ratios at the diameter d., 
corresponding to the input theta locations. The difference between the 
input intensity ratios and those associated with d. is called the residual. 
The program calculates another set of intensity ratios based on 4 new 


diameter d. + Ad., where Ad. is 4 small increment of diameter (usually 
05 microns to start). The intensity ratios calculated using d, + Ad, are 


compared with the input intensity ratios to give a second residual. The 
residuals are compared and program logic determines whether or not the 
Ad increment is producing convergence to the actual diameter. The 
initia] guess diameter is incremented and decremented as necessary until 
a desired level of accuracy is achieved. 

A key to understanding the convergence process is the residual curve 
in Figure 53. It is important that the initial guess, d., be fairly close to 
the actual diameter. 

For example, the residual curve in Figure 33 is for an actual diarneter 
of / microns. If d. is greater than 9 microns, the program will not 
converge to the correct diameter, but rather to a diameter just above 16 
microns. 


It is hard to define exactly the limits of d. for any given actual 


diameter. Figure 541s a residual curve for 5 microns and shows an upper 


oT 


limit of 9microns for d. A thorough study to define limits ford. has 
not been conducted, although convergence has always been attained by 


guessing d, within +1 micron for diameters ranging from 5 to 9 microns. 


RESIDUAL VS DIAMETER 


*10° 
30.0 


RESIDUAL 
10.0 15.0 20.0 25.0 


5.0 


0.0 





0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 {0.0 {1.0 {2.0 
DIAMETER (METERS) ar 


Figure 53. Residual Versus Diameter for 71m Fiber 
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RESIDUAL VS DIAMETER 





6.0 7.0 8.0 9.0 10.0 11.0 12.0 13.0 14.0 15.0 
DIAMETER (METERS) *10° 


Figure 54. Residual Versus Diameter for Sjim Fiber 
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3. EXPOSURE 

EXPOSURE simulates the operation of the Micron Eye Dy taking 4 
photograph of intensity ratio data. EXPOSURE reads in the perfect data 
generated by DATAMAKR and digitizes the data with respect to a threshold 
intensity ratio. 

The optimum intensity ratio for an exposure is termed the 
threshold intensity ratio. The threshold intensity ratio is located 
where the absolute value of the derivative of the intensity ratio with 
respect to theta is a maximum. This assures that the threshold intensity 
ratio is located in a region of the curve which is closest to 4 straight line. 
This reduces the effects of subsequent interpolation errors. 

The threshold intensity ratio is calculated by the subroutine 
DERIY. Since the Input intensity profile curve will have at least three 
maximum derivative points, DERIVY calculates the average of the three 
Intensity ratios. This average is then passed to the calling program, 
EXPOSURE. 

EXPOSURE searches the intensity ratios and compares them with 
the value of the threshold intensity ratio. This process corresponds to the 
Micron Eye addressing 4 pixel and comparing its voltage to the threshold 
voltage. If an intensity ratio is greater than the threshold value, the ratic 
Is assigned a digital value of 1. If the intensity ratio is less than the 


threshold value, the ratio is assigned 4 digital value of 0. 
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‘The digitized data is then searched to find the theta locations at 
which the digital intensities change from 9 --> 1, or from 1--> 0. The 
theta locations are averaged to produce a theta location which is very 
clase to the threshold intensity. Because this averaging process is 
occuring in the regions of steepest slope, error is assumed to be 
minimized since the curve can be approximated as a straight line (refer to 
Figure 30 in Chapter 6). 

EXPOSURE provides the user with the option of introducing error 
into the perfect data. The program prompts the user to input the desired 
error and calls the subroutine RANDOM. RANDOM introduces randorn error 
into the intensity ratios and returns the imperfect data to EXPOSURE. 
EXPOSURE digitizes the imperfect data and searches for the occurences 
where 9-->1 and 1-->0. Imperfect data will have many more occurences 
of the digital intensities changing from 1-->9 and O-->1. EXPOSURE will 
average the theta locations if they are in the same location (i.e., within 
+.005 radians). This method has been tested for random error up to 5 


percent. 
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